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We study the gluon and ghost propagators of lattice Landau gauge in the strong-coupling limit 
j3 — in pure SU(2) lattice gauge theory to find evidence of the conformal infrared behavior of 
these propagators as predicted by a variety of functional continuum methods for asymptotically 
small momenta q 2 <C Aq Cd . In the strong-coupling limit, this same behavior is obtained for the 
larger values of a 2 q 2 (in units of the lattice spacing a), where it is otherwise swamped by the gauge 
field dynamics. Deviations for a 2 q 2 < 1 are well parameterized by a transverse gluon mass oc 1/a. 
Perhaps unexpectedly, these deviations are thus no finite-volume effect but persist in the infinite- 
volume limit. They furthermore depend on the definition of gauge fields on the lattice, while the 
asymptotic conformal behavior does not. We also comment on a misinterpretation of our results by 
Cucchieri and Mendes in Phys. Rev. D81 (2010) 016005. 
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I. INTRODUCTION 

The Green's functions of QCD are the fundamental 
building blocks of hadron phenomenology More- 
over, their infrared behavior is also known to contain es- 
sential information about the realization of confinement 
in the covariant formulation of QCD, in terms of local 
quark and gluon field systems. In relation to the gluon 
and ghost propagators of Landau gauge QCD the Dyson- 
Schwinger equation studies of Refs. |5|,|6| established that 
the gluon propagator alone does not provide long-range 
interactions of a strength sufficient to confine quarks, 
which dismissed a widespread conjecture from the 1970's 
going back to the work of Marciano, Pagels, Mandelstam 
and others. It was concluded that the infrared domi- 
nant correlations are instead mediated by the Faddeev- 
Popov ghosts of this formulation, whose propagator was 
then found to be infrared enhanced, in agreement with 
the Kugo-Ojima confinement criterion and thus consis- 
tent with the conditions for confinement in local quantum 
field theory @,0, |]. 

This infrared behavior was subsequently confirmed by 
a variety of studies based on different functional methods 
in the continuum which all led to the same conformal in- 
frared behavior for the gluonic Green's functions of Lan- 
dau gauge QCD. These include studies of their Dyson- 
Schwinger Equations (DSEs) [8] , Stochastic Quantization 
Q , and of the Functional Renormalization Group Equa- 
tions (FRGEs) [3. In fact, this conformal infrared be- 
havior of gauge-field correlations in non-Abelian gauge 
theories with confinement is directly tied to the validity 
and applicability of the framework of local quantum field 
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theory to such theories. 

However, with the notable exception of pure SU(2) lat- 
tice gauge theory in two dimensions [ll[ , this is not what 
is being observed with current lattice implementations of 
Landau gauge [p3 - [i~9l ]. Rather, these are much more con- 
sistent with an essentially free ghost propagator together 
with a massive and hence non-vanishing gluon propaga- 
tor in the infrared, in qualitative agreement with DSE 
solutions proposed in the studies of Refs. (2fj| - [26l |. Note 
that the functional equations of continuum quantum field 
theory admit both types of solutions, the scaling solution 
with the predicted conformal infrared behavior and the 
decoupling solution with a massive gluon propagator [27l | . 
At present, the appealing aspects for fundamental rea- 
sons of the scaling solution may be seen to stand against 
an overwhelming evidence for the massive one, driven by 
the numerical results from present lattice implementa- 
tions of the Landau gauge. We will review the situation 
briefly in the next section. 

In this paper we study the gluon and ghost propaga- 
tors in the strong couplin g lim it, /3 — > 0, of pure SU{2) 
lattice Landau gauge [28L l29l]. This limit can be inter- 
preted as the limit of infinite lattice spacing, a — ¥ 00, at 
a fixed physical scale as set, e.g., by the string tension 
which behaves as a 2 a oc — ln/3 for j3 — > 0. Alternatively, 
when considering all momenta in lattice units, q ~ 1/a 
as we do here, the same strong-coupling limit can also 
be interpreted as a hypothetic limit in which all physi- 
cal momentum scales such as the string tension, Aqcd or 
the lowest glueball mass are sent to infinity, i.e., formally 
as a — y 00 or Aqcd — ► 00. Both interpretations are of 
course equivalent. Either way, all momenta and masses 
in lattice units 1/a are infinitely small relative to the 
physical scale of the theory which is precisely what we 
need for an analysis of the asymptotic infrared behavior 
of its correlation functions. 

By definition, there is then no scaling limit. All mo- 
menta are zero in physical units. Having said that, we 
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will nevertheless distinguish large lattice momenta, with 
a 2 q 2 >• 1, from intermediate, a 2 q 2 ~ 1, and low mo- 
menta, a 2 q 2 <C 1, in the strong-coupling limit. Physi- 
cally, they are all in the asymptotically-far infrared, but 
for sufficiently large lattices, L/a 1, only small mo- 
mentum modes will be affected by finite-size corrections 
which can be assessed by varying the lattice size, L/a. 

With this in mind, we will show that the conformal 
infrared behavior is indeed observed at large lattice mo- 
menta in the strong-coupling limit. Perhaps surpris- 
ingly, however, the deviations at intermediate and small 
momenta reveal only a negligible (though systematic) 
lattice-size dependence. In agreement with the standard 
lattice Landau-gauge simulations at finite /3, we find that 
the small momentum behavior is well parameterized by 
a transverse gluon mass. This is not predominantly due 
to the finite volume, but it depends on the actual lat- 
tice discretization of the gauge fields, while the large 
momentum-behavior does not. So the latter is consis- 
tent with scaling while the data at small momenta shows 
decoupling, albeit being discretization dependent. 

The intermediate momentum region is characterized 
by the transition between scaling and decoupling around 
a 2 q 2 « 1. This transition sets a scale on an infinite lat- 
tice with a/L — > 0, which corresponds to a transverse 
gluon mass aM of the order one in the strong-coupling 
limit. This scale is thus not related to the string tension 
or glueball masses in the strong-coupling regime which 
all behave as am oc — In /3 for /3 — > 0. The occurance 
of such a new infrared scale aM ~ 1 seems artificial 
and its discretization dependence might in fact signal the 
breakdown of Slavnov- Taylor identities for minimal lat- 
tice Landau gauge in the strong-coupling regime. 

The paper is organized as follows: In Sec. [TT] we re- 
view the expectations for the infrared behavior from the 
continuum studies as background and motivation. The 
different parts of this section may be consulted individu- 
ally or skipped on a first reading as convenient. In Sec. Mil 
we show that our numerical results for the Landau-gauge 
gluon and ghost propagators in the strong-coupling limit 
do in fact show the scaling behavior for lattice momenta 
with a 2 q 2 S> 1, that the critical exponent can be ex- 
tracted in good agreement with continuum predictions 
from this data, and that the deviations from conformal 
scaling for a 2 q 2 < 1 are well parameterized by a trans- 
verse gluon mass M « l/a in the infinite- volume limit. 
In Sec. lIVI we compare various lattice definitions of gauge 
potentials, which would all be equivalent in the contin- 
uum limit, and show that the essential features of the 
scaling branch at large lattice momenta, such as crit- 
ical exponent and coupling, are unique. We further- 
more demonstrate that the massive branch observed for 
a 2 q 2 < 1 does depend on the lattice definition of the 
gluon fields, and that it is thus not unambiguously de- 
fined. We interpret this as an ambiguity in the definition 
of Landau gauge on the lattice which precludes a corre- 
sponding definition of a measure for gauge-orbit space in 
presence of Gribov copies [3(| ■ One might still hope that 



this ambiguity will go away at non-zero /? in the scaling 
limit. While this is true at large momenta, we demon- 
strate in Sec. |V]that the ambiguity is still present in the 
low-momentum region, at least for commonly used val- 
ues of the lattice coupling such as (3 = 2.3 or /3 = 2.5 
in SU(2). Our Summary and conclusions are provided 
in Sec. I VII and further technical details are given in two 
appendices. 



II. INFRARED SCALING VERSUS 
DECOUPLING 

The Landau-gauge gluon propagator, in (Euclidean) 
momentum space, is parameterized by a single dressing 
function Z, 

DfM = ^ («W - *f>j £Gp , (1) 

and the ghost propagator by a corresponding dressing 
function G. 



D a G b (p) = 



ab G(p 2 ) 



V 



(2) 



For their infrared behavior, i.e., that of Z(p 2 ) and G(p 2 ) 
for p 2 — > 0, we consider the two possibilities described in 
the subsections that follow. 



A. Scaling 

The prediction of [BNTolj amounts to infrared asymp- 
totic forms 



Z(p 2 ) ~ (p 2 /A| CD ) 2Kz , 
G( P 2 ) ~ ( P 2 /A 2 QCD )- KG , 



(3a) 
(3b) 



for p 2 — > 0, which are both determined by a unique crit- 
ical infrared exponent 



(4) 



with 0.5 < k < 1. Under a mild regularity assumption 
on the ghost-gluon vertex joL the value of this exponent 
is furthermore obtained as [H, Q 



= (93 - Vl201)/98 fa 0.595 



(5) 



The conformal nature of this infrared behavior in the 
pure Yang-Mills sector of Landau gauge QCD is evi- 
dent in the generalization to arbitrary gluonic correla- 
tions 31|: a uniform infrared limit of one- particle irre- 
ducible vertex functions r m ' n with m external gluon legs 
and n pairs of ghost/anti-ghost legs of the form 



(p 2 /Al~- i("- m ) K 



l QCDJ 



(6) 
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when all pf oc p 2 — > 0, z = 1, . . . 2n + to. In particular, 
the ghost-gluon vertex is then infrared finite (with n = 
m = 1) as it must [32j , and the nonperturbative running 
coupling introduced in Q via the definition 



a s (p 2 ) = »-Z{p 2 )G 2 {p 2 ) 



(7) 



approaches an infrared fixed-point, a s — > a c for p 2 — > 0. 
If the ghost-gluon vertex is regular at p 2 =0, its value is 



8tt T 2 (k - l)r(4- 2k) 
N~ c T 2 (-k)T{2k- 1) 



Nr. 



x 0.99 



(8) 



Comparing the infrared scaling behavior of DSE and 
FRGE solutions of the form of Eqs. ([3]), it has in fact 
been shown that in presence of a single scale, the QCD 
scale Aqcd > the solution with the infrared behavior (BT) 
and (J6j) , with a positive exponent K, is unique (33l . l34| 
and nowadays being called the scaling solution. 



B. Decoupling 

This uniqueness proof does not rule out, however, the 
possibility of a solution with an infrared-finite gluon 
propagator, as arising from a transverse gluon mass M, 
which then leads to an essentially free ghost propagator, 
with the free massless-particle singularity at p 2 = 0, i.e., 



Z{p 2 ) ~ p 2 /M 2 , and GO 2 ) ~ const. 



(9) 



for p 2 — > 0. The constant contribution to the zero-mo- 
mentum gluon propagator, D(0) — 1/M 2 , thereby neces- 
sarily leads to an infrared constant ghost renormalization 
function G. This solution corresponds to Kz = 1/2 and 
kq = 0. It does not satisfy the scaling relations (j4]) or 
©. This is because in this case the transverse gluons 
decouple for momenta p 2 <C M 2 , below the independent 
second scale given by their mass M. It is thus not within 
the class of scaling solutions considered above, and it is 
termed the decoupling solution in contradistinction. The 
interpretation of the renormalization group invariant (7) 
as a running coupling does not make sense in the infrared 
in this case, in which there is no infrared fixed-point and 
no conformal infrared behavior. 



C. Continuum versus Lattice Studies 

The functional equations of continuum quantum field 
theory admit two types of solutions. With the current 
implementations of Landau gauge on the lattice strong 
support is provided for the massive decoupling solution 
from simulations. Is this the end of the story? Superfi- 
cially one might say that lattice gauge theory provides an 
ab initio framework whereas functional continuum meth- 
ods yield ambiguous results. We would like to argue that 



this conclusion might be a bit premature, however. First, 
lattice simulations must necessarily be done in a finite 
volume where, strictly speaking, a conformal behavior 
cannot be observed, certainly not for the lowest momen- 
tum values. What is necessary to observe an at least 
approximate conformal behavior of the correlation func- 
tions in a finite volume of extend L, is a wide separation 
of scales, 



n/L < p < Aqcd 



(10) 



such that a reasonably large number of modes with mo- 
menta p sufficiently far below the QCD scale Aqcd are 
accessible 1 whose corresponding wavelengths are at the 
same time much shorter than the finite size L. 2 Sec- 
ondly, gauge-fixing in presence of Gribov copies is not so 
well understood on the lattice either and the relation be- 
tween present implementations of lattice Landau gauge 
and covariant gauge theory in the continuum with lo- 
cal Becchi-Rouet-Stora-Tyutin (BRST) symmetry is far 



from clear after all [37|. Last but not least, the func- 



tional continuum methods are not as ambiguous as they 
might at first appear. At least technically the fact they 
admit both kinds of solutions is not surprising because 
it is well known that these methods have to be supple- 
mented by additional boundary conditions [8j, [9j, l27| . It 
is these boundary conditions that determine which of the 
two solutions will be obtained, and it is essentially these 
boundary conditions in which the two classes of contin- 
uum studies mentioned above differ from one another. 
In particular, the scaling solution requires the boundary 
condition for the subtraction in the ghost DSE such that 



G _1 (p 2 ) -> 0, forp 2 -> 0, 



(11) 



which in Landau gauge then implements both, the unbro- 
ken global gauge charges of the Kugo-Ojima confinement 
criterion and the horizon condition of the original Gribov- 
Zwanziger framework [30l [38[ , by the infrared dominance 
of ghosts. This implies that the Kugo-Ojima confinement 
criterion cannot be derived from DSEs alone. It is implic- 
itly implemented by the boundary condition which 



1 The relevant scale here is that of the minimal MOM (MM) 
scheme 135]. Its relation to the MS scheme, e.g., for zero fla- 
vors and N c = 3, with A-^ 240 MeV yields A MM s« 450 MeV. 
The enormous challenge for lattice simulations thus is to satisfy 
tt/L <C p <C 450 MeV in the scaling region. 

2 Finite-volume effects have been analyzed carefully for the DSE 
scaling solution on a 4-dimensional torus which showed a quite 
significant volume dependence |36| . In particular, it was con- 
cluded that volumes of about 15 fm in length are needed before 
even an onset of the leading infrared behavior can be observed, 
and that up to 40 fm might be required for a reliable quantitative 
determination of an approximately conformal infrared behavior 
from a suitable range of momenta satisfying dlOII . That this de- 
pendence does not match up with that on a lattice became clear 
at the 2007 Lattice Conference, where several results were re- 
ported from large lattices [1341611 , the largest for SU (2) with up 
to 27 fm in size [15], showing practically no tendency to follow 
the predicted finite-size corrections to the DSE scaling solution. 
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leads to the conformal infrared behavior | : 8| , while other 
choices genuinely lead to the massive one [271 ]. 

More generally, from the functional equations alone 
(whether DSEs or FRGEs or both together) both solu- 
tions are possible, in principle. Because only the scaling 
solution is consistent with the conditions for confinement 
in local quantum field theory, based on the cohomology 
construction of a physical Hilbert space over the indefi- 
nite metric spaces of covariant gauge theory from BRST 
symmetry, this appears to be the physically relevant so- 
lution within this framework and it therefore received a 
lot of attention in the functional continuum studies. 

The decoupling solution ([5]) has received renewed at- 
tention [22j-|26| mainly because this is what is being ob- 
served in Landau gauge implementations on the lattice 
[T2l - [l9j . The numerical procedures on the lattice are 
thereby based on minimizations of a gauge-fixing po- 
tential with respect to gauge transformations. To find 
absolute minima is not feasible on large lattices as this 
is a non-polynomially hard computational problem. One 
therefore settles for local minima which in one way or an- 
other, depending on the algorithm, samples gauge copies 
of the first Gribov region. This is mimicked in the con- 
tinuum by the inclusion of Zwanziger's horizon func- 
tional to suppress the gauge copies outside the first Gri- 
bov region within the Gribov-Zwanziger framework. To 
make this framework compatible with the decoupling so- 
lution ^ one then introduces an additional mass term 
[22|,[39j]. While the renormalizability is maintained, one 
no-longer has an exact local BRST symmetry in this 
framework which leads to unitarity violations when at- 
tempting a BRST cohomology construction of a phys- 
ical Hilbert space. This so-called soft BRST breaking 
only matters in the non-perturbative regime and it re- 
lates to the sampling of Gribov copies. A similar effect, 
in terms of reweighting copies inside and outside the first 
Gribov horizon, is achieved by introducing an explicit 
Curci- Ferrari mass [40]. In this case the reweighting is 
explicitly controlled by the mass and results in an ex- 
plicit BRST breaking proportional to that mass. It ap- 
pears that a reweighting of Gribov copies generally causes 
BRST breaking. At the moment it seems questionable to 
us whether one can have a gluon mass without reweight- 
ing of Gribov copies and BRST breaking. 

Within the framework of local quantum field theory, 
which however requires an unbroken BRST symmetry 
with nilpotent BRST charge, the decoupling solution © 
is realized if and only if it comes along with the Higgs 
mechanism. The Kugo-Ojima confinement criterion and 
the infrared scaling of Landau gauge Green's functions 
as a consequence of this criterion cannot be dismissed 
from lattice simulations without a proper definition of a 
non-perturbative BRST symmetry. It is worth remem- 
bering, however, that the apparent ambiguity arises only 
when comparing the gauge dependent Green's functions 
of either approach. By construction, physical observables 
remain of course unaffected by this problem with BRST 
in minimal lattice Landau gauge. One example is the 



Polyakov-loop potential of the pure gauge theory whose 
center symmetry can be used to define an alternative 
confinement criterion which is in fact satisfied by the de- 
coupling solution as well [4l| , regardless of the realization 
of BRST symmetry. 

The agreement between lattice Landau gauge and con- 
tinuum (decoupling) results, i.e., when the restriction to 
the first Gribov region is implemented, seems quite con- 
vincing. The fact that the accuracy and conclusiveness 
of lattice results together with our understanding of the 
functional methods based on local quantum field theory 
have unveiled the conflict between the observed dynam- 
ical gluon mass and BRST symmetry is a great achieve- 
ment. We believe that it will eventually allow us to un- 
derstand the relation between Gribov copies and BRST 
symmetry. 

Meanwhile, the strong-coupling limit of lattice Landau 
gauge provides a powerful tool to study this ambiguity. 
In this unphysical limit, the gluon and ghost correlations 
are solely driven by what should correspond to a non- 
perturbative measure of gauge-orbit space in the Landau 
gauge, i.e., the gauge-fixing and Faddeev-Popov parts of 
the lattice measure. It is this measure that is being as- 
sessed when the gauge-field dynamics is switched off. 

Ghost dominance, the essential condition for the con- 
formal infrared behavior, is then implemented by hand, 
and if there is such a behavior, it should be seen at least in 
this limit in which all momenta are asymptotically small 
(in physical units). Because the strong-coupling limit can 
be interpreted as the formal limit Aqcd — > oo, it is par- 
ticularly well suited to assess whether the conformal be- 
havior of the scaling solution is seen for the larger lattice 
momenta, with a 2 q 2 ir 2 a 2 /L 2 and thus well clear of 
finite-size effects, after the upper bound in ([TU| has been 
removed. This range of lattice momenta would otherwise 
be dominated by the dynamics due to the gauge action 
whose presence thus obscures any potentially conformal 
behavior there. As we demonstrate below, its removal 
does in fact reveal a scaling behavior ([3]), (Q} for the first 
time in lattice simulations, and it furthermore allows to 
study the corrections to scaling from finite- volume effects 
on reasonably small lattices with some systematics. 

It follows unambiguously from this study that finite- 
volume effects play a minor role, however, and that up to 
these small effects, decoupling is observed at small lattice 
momenta a 2 q 2 <C 1, with a mass parameter M cx 1/a. 
But the strong-coupling limit also serves to isolate a 
discretization ambiguity which manifests itself in depen- 
dences on the lattice definition of gauge fields underlying 
the respective lattice Landau gauges and their measures. 
While this noticeably affects the decoupling branch at 
a 2 q 2 < 1, the critical exponent and coupling of the scal- 
ing branch at large a 2 q 2 are rather insensitive to the dis- 
cretization. 3 



3 Note, however, that the extraction of a scaling exponent from the 
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FIG. 1: The gluon propagator (left) and its dressing function Z (right) versus lattice momentum x = a 2 q 2 for different lattice 
sizes in the strong-coupling limit. For illustration purposes the data for Z (right) is compared to the continuum predictions 
from decoupling (solid) and scaling (dashed) in infinite- volume limit. That is, the respective exponents have not been fitted 
but set to the expected decoupling and scaling values k = 0.5 and k ~ 0.595. 



III. STRONG-COUPLING LIMIT OF 
STANDARD LATTICE LANDAU GAUGE 



stopping criterion 



We simulate pure 5/7(2) gauge theory in the strong- 
coupling limit by generating random link configurations 
{U}. These are sets of 5/7(2) gauge links, 



XLL^ 



(12) 



equally distributed over (it , u) x ^ G S 3 . Those configura- 
tions are then fixed to the standard lattice Landau gauge 
(SLG) using an over-relaxation algorithm that iteratively 
minimizes the SU (2) gauge-fixing functional, for SLG, 



e := max tr [(V^ M )(V^)] < 10" 



13 



(15) 



is satisfied for every site on the lattice. Gluon and ghost 
propagators are then calculated in momentum space em- 
ploying standard techniques. The dressing functions, Z 
and G, are extracted from the known tree-level form 
[Eqs. (|A6|) and (|A11[) ] of the respective lattice propa- 
gator, see Appendix lAl for further details. 



Vu[g} = A^2(l 



(13) 



X.fi 



where the = gxU x ^g\^^ are the gauge-transformed 
links. The Landau gauge condition for the stationar- 
ity of Vj/[<7] under gauge transformations g is given by 
the lattice divergence, F X (A 9 ) = V*Ag M = 0, where 
denotes the lattice backward derivative and A 9 xfl is the 
lattice gluon field of SLG, defined by 



(14) 



in terms of the gauge-transformed link U 9 . To imple- 
ment the minimal Landau gauge with a sufficient accu- 
racy the over-relaxation algorithm is iterated until the 



0.64 



0.60 



0.56 




L/a 



ghost propagator is affected b y th e Gribov ambiguity as shown 
in the follow-up study of Ref. [42j| . 



FIG. 2: The zero-momentum limit (|17l) of the strong-coupling 
gluon propagator over L/a from global fits of the form (|18c[) . 
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FIG. 3: x 2 / n df an d K z values of fits to the gluon propagator data for different lattice sizes and different lower bounds xi. The 
upper bound x u = 14 has been kept fixed. The respective fit functions are given in the upper panels. 



A. Gluon propagator 

In Fig. Q] the data for the gluon propagator, 

D(x) = Z(x)/x , with x = a 2 q 2 , (16) 

and its dressing function Z are plotted against the lattice 
momenta a 2 q 2 denned in Eq. (|A7[) . The propagator is 
observed to increase with momentum, while it plateaus 
at low momenta. Perhaps unexpectedly, however, this 
happens irrespective of the lattice size (N — L/a) at 
around a 2 q 2 rs 1. It is therefore not primarily due to the 
finite volume. Rather, on sufficiently large lattices the 
observed mass behaves as 

M 2 = limD -1 (a;) oc l/a 2 (17) 

x—>-0 

in the strong-coupling limit with hardly any significant 
dependence on L. In particular, if there is a systematic L 
dependence at all, the zero momentum limit of the gluon 
propagator tends to slowly increase with the volume as 
shown in Fig. [5] It certainly extrapolates to a finite value 
oc 1 1 a 2 in the limit L/a — > oo . 

In the right panel of Fig. [T] we furthermore compare the 
strong-coupling data for the gluon dressing function Z, 
to the predicted forms corresponding to decoupling ©, 
with Zd(x) = CdX, and scaling (I3al) . with Z s [x) = c s x 2kz 
where the value of kz is not fitted but taken from Eqs. (HJ 
and ©, = 0.595, for comparison. 

With Cd and c s appropriately adjusted to the 56 4 data, 
we find that the decoupling solution provides a very good 
description of the low momentum region, while the large 
momentum branch approaches the scaling solution with 



an exponent kz clearly above 0.5 (the fitting procedure 
described below leads to a conservative estimate of about 
kz = 0.57(1) for the 56 4 lattice, for example). 

In order to assess the asymptotic form at large lattice 
momenta x = a 2 q 2 more quantitatively, we have fitted 
the gluon propagator data to the following three forms: 

D a (x) =cx 2Kz ~ 1 , (18a) 

D b {x) =cx 2Kz ~ 1 + d, (18b) 

D c {x) = c{d + xf Kz - 1 . (18c) 

D a describes pure scaling with an effective exponent 
while Di, and D c accommodate the transition between 
decoupling at small x and scaling at large x in different 
ways. In order to analyze the scaling exponent kz we 
have used all three forms to fit the data in various win- 
dows with increasing lower bound xi. The results are 
fairly insensitive to variations of the upper bound x u in 
some range sufficiently close to the maximum value of 
x = 16. We used x u = 14 in all fits. 

The results of these fits for kz as functions of the lower 
bound xi with the corresponding x 2 /ndf are summarized 
in Fig. [3] The pure scaling model D a (x) in (I18ap cannot 
describe the full momentum range but leads to good and 
stable fits to the data for xi > 1, the form D c {x) in (|18c[) 
provides the best global description of the data over the 
full momentum range. For xi between 1 and 3 the results 
for kz from all three models are consistent with each 
other within errors and all with x 2 /ndf of around 1. The 
values of kz that result for the different lattice sizes from 
all three fit models with x; = 1 are shown in Fig. |U 

As before, these values show very little systematic de- 
pendence on L/a. A slight tendency to drift towards 
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0.65 



0.60 



0.55 - 



0.50 - 



D(x) = c(d + x) 1 



12 



24 



32 
L/a 



56 



FIG. 4: The infrared exponent kz for lattice sizes L/a = 
12, 24, 32 and 56 as obtained from the strong-coupling gluon 
data on x € [1, 14] with the three different fit models in (|18l) . 



larger values in larger volumes is observed, but this might 
not be a significant effect. Assuming that there is no lat- 
tice size dependence to fit the values from model D c (x) in 
(I18c[) for the 4 different lattice sizes by a constant yields 
an average of kz — 0.563(1), consistent with a global av- 
erage over all values in Fig. |4j while the 56 4 data alone 
with the same model gives kz = 0.564(1). For compar- 
ison, model Dt,(x) in (|18bp . with the largest errors, for 
the 56 4 data on x € [1, 14] yields k z = 0.568(9). 



B. Ghost dressing function 

We performed similar fits to extract the exponent kg 
from the strong-coupling ghost dressing function G as 
shown in Fig. [5] Those fits are less robust with a more 
pronounced systematic uncertainty due to the fit model 
dependence. This is mainly because of the wider transi- 
tion region, from G = const, at small x to G ~ x~ K ° at 
large x, which is under less control here. In fact, a pure 
power law is at best observed only for the very largest val- 
ues of the lattice momentum, in the range above x w 10 
or so. We again used three different fit models analogous 
to those for the gluon data in (JTHJ) , one describing pure 
scaling at large lattice momenta and two that interpo- 
late between decoupling © and scaling (l3b|) . This time, 
however, we fit the inverse of the ghost dressing function 
G as follows: 



G~ 1 (x)=cx KG , 
G^ l {x) = cx KG +d, 
G~ l {x) =c{d + x) KG 



(19a) 
(19b) 
(19c) 



With the same method as used above, the pure scaling 
form G a with lower bounds xi around 12 leads to values 
of kg around 0.52 for the 56 4 lattice with a tendency to 
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FIG. 5: The ghost dressing function G for different lattice 
sizes in the strong-coupling limit compared to the continuum 
predictions from decoupling (solid) and scaling (dashed) with 
the exponent from Eq. (JS| in infinite- volume limit (not fitted). 
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FIG. 6: k versus a/L for the ghost and gluon propagators. 
Grey-colored bands mark the variation of k with the fit model 
and fitting window for either propagator. 



further increase with xi, whereby the \ 2 / n ^ decrease. 
But the quality of those fits is still rather poor indicating 
that an asymptotic scaling cannot be isolated from the 
transition region in the ghost renormalization function. 

The fit models (|19b[) and (|19c[) take this transition re- 
gion better into account. The values for kg from Gb in 
(|19bj) are the most stable ones but tend to be rather large 
with kg around 0.68. Model G c leads to exponents kg 
most consistent with the scaling relation in Eq. (fj|, with 
kq in a range between 0.55 and 0.62 depending on xi. 

The results for both exponents, kz and kg, are sum- 
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marized in Fig. [51 For the gluon exponent Kz these are 
the same data points as in Fig. 0] with the dark grey 
band indicating the errors and systematic uncertainties 
due to the fit model. The kq values are shown for a 
variety of lower bounds xi as just described with model 
(|19cp overlapping the gluon results. The light grey band 
in Fig. [6] indicates the considerably larger uncertainties 
in the ghost exponent, which are mainly the systematic 
ones due to the particular difficulty in modelling the wide 
transition region between decoupling and scaling there. 

Nevertheless, the set of all values extracted for kg from 
the 56 4 data center around kg = 0.60(8) which includes 
the range for the gluon exponent and is thus fully con- 
sistent with the scaling relation kz = kg in Eq. (g]). 



IV. DIFFERENT GAUGE-FIELD DEFINITIONS 
ON THE LATTICE 

Strong-coupling configurations are very rough with 
links distributed uniformly over the parameter space, the 
3-sphere for SU(2). The strong-coupling limit is there- 
fore an ideal testbed for different lattice definitions of 
gauge-fields which correspond to different choices of co- 
ordinates that agree only near the identity, or in the con- 
tinuum limit. The standard definition (|14[) for exam- 
ple corresponds to choosing separate coordinates for the 
Northern (NH) and Southern Hemispheres (SH) of S 3 
in the case of SU(2). Strictly speaking, the SLG gluon 
propagator therefore corresponds to an average for each 
link of the contributions from NH and SH to the expec- 
tation value in Eq. (|A5I) . 

The maximal chart is provided by stcrcographic pro- 
jection which covers the whole sphere except for the 
South Pole. A definition of SU (2) gauge fields on the 
lattice based on stereographic projection is possible as 
follows, 



where 



i 



2 ^ xfi 



u 



X f t 



(20) 



(21) 



It agrees with the standard definition near the North 
Pole, and in the continuum limit, but the South Pole 
is now at infinity and the gauge fields in (|2"0)) are thus 
non-compact variables. The associated Landau g aug e is 
the modified lattice Landau gauge (MLG) of Ref. [4J|. It 
follows from the stationarity condition of the modified 
gauge-fixing functional, 



V u [g]=-8j2 ln (^ + i trU ^j 



1 , 1 



(22) 



X , ft 



with respect to gauge transformations g 1 and reads 

Fx{A 9 ) = F(A?) = V»l*„ = . (23) 



When comparing MLG to the ever popular SLG, there is 
no advantage that the SLG has over the MLG. A promis- 
ing particular feature of the MLG on the other hand is 
that it provides a way to perform gauge-fixed MC simu- 
lations sampling all Gribov copies of either sign (of the 
Faddeev-Popov determinant) in the spirit of BRST. This 
feature will be explored in a forthcoming study. Here we 
simply use the MLG for comparison in the standard way, 
i.e., we gauge- fix configurations via minimization of the 
MLG functional in Eq. (JUJ). The F-P operator of the 
MLG is given in (|B7|) in Appendix [B] 

Both lattice definitions of Landau gauge have the same 
continuum limit, and any differences between MLG and 
SLG data at finite lattice spacings are lattice artifacts. 
It is also worth mentioning that gauge configurations 
fixed to MLG do not satisfy the gauge condition of SLG 
and vice versa. Nevertheless, exact transversality, i.e., 
q l _ l (k)A tl (k) — or q^(k)A^{k) = 0, is satisfied at finite 
lattice spacing a for both of them equally, with their re- 
spective Lorenz conditions V b ^A xtl = or V^A xfl = 
and midpoint definition, if the momenta q^ik) are de- 
fined as aq^{k) = 2 sin(7rfc M /A M ) with integer valued 
k,j, £ (—N^/2, Np/2], The identification of this so de- 
fined q^{k) with physical momentum is the usual tree- 
level correction for the Wilson gauge action. It is a spe- 
cial feature of MLG and SLG that their lattice Landau 
gauge conditions define gluon fields that are transverse 
in this physical momentum at any finite lattice spacing. 



A. Gluon propagator 

The data for the gluon propagator of SLG (red filled 
diamonds) is compared to that of MLG (blue filled cir- 
cles) in Fig. [JJ There we also show data for the gluon 
propagator where either 



a A x d l i — Ux^x^" 



(no sum n) 



or 



aA'l 



b»o- a /2 from U x 



(24) 



(25) 



were used to define lattice gluon fields based on the ad- 
joint representation, A adl (black open diamonds), and 
thus blind to the center 44]; or on the tangent space at 
the identity A ]n (green crosses). In these two cases, A ad i 
and A ]n , for the purpose of a qualitative comparison, we 
simply use the gauge configurations of the SLG to calcu- 
late the gluon propagator. Especially for A ln this implies, 
however, that the condition q^(k)A^(k) — is satisfied 
at best approximately and nowhere near the precision of 
the other two (SLG and MLG). The residual uncertainty 
due to other possible tensor structures in the gluon prop- 
agator (|A5|) then causes the somewhat larger errors for 
this definition as seen in Figs. [7] and [S] 

In Fig.jJJwe have first fitted the data from all four def- 
initions to D c in Eq. (I18cp which provides the best over- 
all description in the full momentum range as mentioned 
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FIG. 7: The strong-coupling gluon propagator over a 2 q 2 for 
the various definitions of gauge fields in Q27p . All on 32 4 
lattices and normalized to the scaling branch after fitting to 
D c in (fTgc) ; all with k z = 0.562 from the fit to the SLG data. 



FIG. 8: Same data as in Fig. [7j but fitted to D b in (|18b|) 
with fixed kz = 0.562 and mass term subtracted in order to 
demonstrate how the data from all definitions then collapse 
onto a unique scaling curve a x 2kz ~ 1 . 



above. In order to demonstrate how the other definitions 
compare to the SLG, we keep its value for the exponent 
fixed when fitting the other data, i.e., Kz — 0.562 as 
obtained for L/a = 32 in SLG is used in all fits. Rela- 
tive to the scaling branch oc x 2kz ^ 1 for large x = a 2 q 2 
we then observe a strong definition dependence in the 
(transverse) gluon mass term at small x. The relative 
weight of the two asymptotic branches, scaling at large x 
and massive at small, is clearly discretization dependent, 
and this dependence cannot be compensated by finite 
renormalizations as is manifest in the data of Fig. [7] 

The fact that the observed mass from the zero- 
momentum limit (|17[) behaves as M oc l/a is a first indi- 
cation that it is indeed this massive branch which is the 
ambiguous one. This is consistent with the fact that the 
definitions of gauge fields on the lattice, which agree at 
leading order, all differ at order a 2 . In particular, with 

U = cos ^ + i3 ■ (f sin ^ and A = tr a A , (26) 

for SU(2), the four different gauge field definitions (lf~4l) . 
flU), ([M]) and ([25]) simply correspond to 



a A = 


^2sin| 




(27a) 


a A = 


j) 4 tan j 




(27b) 


aA ad > = 


j> sin^ 




(27c) 


aA 1 " = 


•j) , with <j) — (f>4 




(27d) 



which clearly all agree only at leading order in the limit 
a — > 0. From the order a 2 differences, the corresponding 
Jacobian factors lead to likewise different lattice mass 
counter-terms for each of the 4 definitions. This is well 
known from lattice perturbation theory where the lattice 



Slavnov- Taylor identities guarantee, however, that the 
gluon remains massless at every order by cancellation of 
all quadratically divergent contributions to its self-energy 
for each of the definitions. In the strong-coupling limit of 
minimal lattice Landau gauge, with an effective decou- 
pling mass © behaving as M 2 oc l/a 2 , such a contri- 
bution survives. This contribution depends on the mea- 
sure for gauge fields whose definition from minimal lat- 
tice Landau gauge beyond perturbation theory is there- 
fore ambiguous. The observation that these differences 
matter here explicitly demonstrates the breakdown of the 
lattice Slavnov- Taylor identities in minimal lattice Lan- 
dau gauge in the non-perturbative domain. 

To assess whether this ambiguity has an influence on 
the exponent kz, we have also used the fits to the form 
Df, in Eq. (|f 8b[) . again with kz = 0.562 fixed from the 
32 4 SLG data. This fit model leads to somewhat larger 
X 2 /ndf arising from the transition region around a 2 g 2 ~ I 
which this form does not describe quite as well as D c , see 
the discussion in the previous section. Having obtained 
the fits to Db allows us to subtract the constant d, how- 
ever, which then makes the normalized data of all four 
definitions nicely collapse onto a unique curve oc x 2kz ~ 1 
as seen in Fig. [HI 

The scaling region in the strong-coupling data for the 
gluon propagator from all 4 definitions is fully consistent 
with a unique exponent Kz of around the SLG value with 
a conservative estimate of an infinite- volume extrapola- 
tion of kz = 0.57(3). 



B. Ghost dressing function 

Equally consistent definitions of transverse gauge- 
fields, gauge conditions and Faddeev-Popov operators are 
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FIG. 9: Inverse ghost dressing functions in the strong- 
coupling limit of minimal lattice Landau gauge using standard 
(SLG) and modified (MLG) gauge fields/conditions. 



available for SLG and MLG which allow us to compare 
their respective ghost propagators also. In Fig.|H]we show 
the strong-coupling data for the (inverse) ghost dressing 
functions from these two alternative definitions of lattice 
Landau gauge, again normalized to the scaling branch at 
large momenta. For this normalization we have fitted the 
data to model G c in (119c[) assuming the scaling relation 
kg = K z with the value of kz = 0.562 from the 32 4 SLG 
gluon data. We then again observe that the two defini- 
tions approach each other in the scaling branch oc x~ K 
at large x but deviate in the decoupling branch at small. 
Within the fit model uncertainties the ghost data from 
both definitions is again consistent with a unique scal- 
ing exponent, and with the scaling relation kq = Kz- 
And again the deviations in relative strengths of the two 
branches cannot be compensated by renormalization. 

Note that in both cases, for the strong-coupling gluon 
and ghost propagators, the normalization constants de- 
termined from their respective scaling branches are actu- 
ally not arbitrary. Their product is related to the critical 
coupling © and thus not independent but also unique as 
we will discuss next. This is consistent with the conclu- 
sion that it is the decoupling behavior at low momenta 
which is ambiguous but not the scaling behavior at large. 



V. RUNNING COUPLING 

The renormalization-group invariant product of Lan- 
dau gauge gluon and ghost dressing functions © defines 
a running coupling. Its perturbative behavior can be de- 
termined from unrenormalized bare lattice data for the 
minimal Landau gauge propagators [45 47] . Its relation 
to the running coupling in the MS scheme is known to 
four loops and it can provide a valuable alternative to 



the MS coupling in phenomenological applications (35j j . 
This furthermore permits an independent additional lat- 
tice determination of the QCD scale parameter from 
continuum extrapolation of the bare product of the lat- 
tice Landau gauge propagators |45l - l47| . 

Without renormalization the bare lattice propagators 
are normalized so as to reproduce their respective tree- 
level forms [Eqs. (|A6I) and (jAllI) ] for the trivial link con- 
figuration, when all links U are set to the identity el- 
ement. Then, however, there is no significance in the 
constant prefactors of the individual propagators. This 
is why we removed these overall constants when compar- 
ing SLG and MLG data for these propagators individu- 
ally as in Figs. [7] and O above. Their product (|7J does 
not get renormalized, however, and should therefore be 
independent of the lattice definition used for the gauge 
fields up to discretization errors. We will discuss this 
definition independence separately in the strong-coupling 
limit, where the discretization effects are largest, and at 
finite lattice couplings f3 where it provides an important 
consistency check for the lattice determinations of the 
QCD scale parameter from the large momentum data of 
this strong running coupling as described above and in 
more detail in Ref. 



47] 



A. Strong-coupling limit 

The predicted infrared scaling (J3J with kz = kg im- 
mediately implies that the running coupling defined in 
Eq. approaches an infrared fixed point, a s — > a c for 
p 2 — > 0. Standard continuum conventions of course need 
rescaling g 2 Z — > Z when comparing to lattice defini- 
tions such as (IT4|) . The predicted conformal scaling in 
the strong-coupling limit, with 

Z = c z (a 2 q 2 ) 2fi and G" 1 = c G (a 2 q 2 ) K , (28) 

would therefore imply that the coupling ([7]) should be 
constant with 



(29) 



Note that its value is thus determined precisely by those 
multiplicative constants in the propagators that were ir- 
relevant to the analysis in the previous section. They 
have to be extracted from the bare lattice data without 
rescaling or renormalization. Besides the critical scaling 
exponent n, the critical coupling a c is an independent 
additional prediction from infrared scaling and it is de- 
termined by these constants. 

In complete agreement with the general observation 
of conformal scaling at large momenta in the strong- 
coupling limit, the product Q of the gluon and ghost 
dressing functions levels at a c w 4 for large a 2 q 2 , as 
seen in Fig. [TU] This is just below the upper bound 
a max ra 4 46 for su(2). The fact that a c obtained here 
should be slightly smaller than this maximum value also 
complies with the continuum prediction in Ref. @ : Con- 
sistent with the results of the previous section this is the 
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FIG. 10: The running coupling, a s , for the standard (SLG) 
and modified (MLG) lattice Landau gauge at /3 = 0. The 
dotted line, a™ ax , is the critical coupling in Eq. ([§} for iV c = 2. 



FIG. 11: Data on a a for the standard (SLG) and modified 
(MLG) lattice Landau gauge at /3 = 2.3 using a 56 4 lattice. 
The lines are spline interpolations to guide the eye. 



expected trend for an exponent k which is slightly smaller 
than the value in ([5]). 

It is furthermore quite compelling that this result, a c 
close to 4, is nearly independent of the gauge-field defi- 
nition, likewise. It is almost identical for SLG and MLG, 
see Fig. [TO] Predominantly driven by the ghost propa- 
gator, the violations to conformal infrared scaling in the 
form of a momentum dependence of a s in the strong- 
coupling limit, set in as soon as the ambiguity in the 
definition of minimal lattice Landau gauge does. 



B. Intermediate lattice couplings 

The significant differences observed at small momenta 
between SLG and MLG so far were linked to the strong- 
coupling limit in which discretization effects are enhanced 
to the extreme. Even though these effects might be ex- 
pected to disappear in the continuum limit, eventually, 
it is important to assess to what extent they survive at 
finite /3. This is of relevance especially to the determi- 
nations of the QCD scale parameter A^g based on this 
running cou pling in the perturbative domain from lattice 



simulations [45H47| . In particular, it is of paramount im- 



portance to verify that this ambiguity of minimal lattice 
Landau-gauge vanishes there, or else to find ways to in- 
clude it in the estimate of the systematic uncertainties, 
if necessary. 

As a first check, we have performed simulations in 
SU(2) where gauge configurations were generated with 
the one-plaquette Wilson action at (3 = 2.3 and 2.5. The 
configurations were then fixed, as above, to SLG and 
to MLG, respectively, to measure and compare a s on 
those two sets. For (3 — 2.3 the results are shown in 
Fig. 1111 Luckily, for the a s project mentioned above, we 
find no significant deviations between the two definitions 



at large momenta. Both agree within the statistical er- 
rors there, and the ambiguity does not appear to affect 
the high momentum behavior of a s (q 2 ) as relevant to 
the fits of its perturbative expansion to extract the QCD 
scale. Of course this should nevertheless be examined 
more carefully as one of the possible systematic uncer- 
tainties also for the relevant 5/7(3) configurations with 
dynamical quarks in future studies. 

At lower momenta one observes significant differences 
between SLG and MLG, especially in the transition re- 
gion around the maxima at a 2 q 2 rs 0.25. In units of the 
string-tension a, with a 2 a = 0.145 for j3 — 2.3 from [ic| . 
this transition region corresponds to q 2 ~ 1.7a, or with 
a = (440 MeV) 2 , in physical units, q « 570 MeV at max- 
imum. These observed differences persist for /3 = 2.5. A 
careful comparison of the strength of the effect for dif- 
ferent /3 but comparable physical volumes, to disentangle 
discretization and finite-volume effects at low momenta, 
is left for future studies, however. 

Because the Gribov copies of SLG and MLG differ, it 
seems quite plausible that the observed ambiguity in the 
minimal lattice Landau gauge beyond perturbation the- 
ory is closely related to the Gribov-copy problem. This 
is also suggested by our data: the effect is predominately 
due to the ghost propagator which is known to be affected 
by this problem at low momenta [42], l49l - [5ll | . 



VI. SUMMARY AND OUTLOOK 

We have studied gluon and ghost propagators of 
pure SU (2) minimal lattice Landau gauge theory in the 
strong-coupling limit. This unphysical limit probes the 
gauge field measure of the minimal lattice Landau gauge 
for there is no contribution from the Yang-Mills (pla- 
quette) action. The Faddeev-Popov determinant is im- 
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plicitly included by collapsing the gauge orbits onto the 
first Gribov region as sampled by the minimal Landau 
gauge implementation used on the lattice. The strong- 
coupling limit can therefore be thought of as a means to 
implement, by hand, a lattice analogue of the infrared 
dominance of ghost contributions in functional methods 
such as DSE or FRGE studies. 

As expected in the formal limit Aqcd — > oc, it is then 
observed that the dressing functions of both propagators, 
Z and G, show the conformal scaling behavior 

Z oc (aV) 2K and G oc {a 2 q 2 )- K 

for large lattice momenta, a 2 q 2 1, well clear of the re- 
gion where finite-size effects should be expected. These 
effects, on the other hand, turn out to be surprisingly 
small and the combined gluon and ghost data is consis- 
tent with an L/a -> oo extrapolation of a critical ex- 
ponent K — 0.57(3). This scaling branch at large a 2 q 2 
furthermore leads to a critical coupling of a c ~ 4 which 
is just below the predicted maximum a™ ax s» 4.46 for 
SU(2). These results show very little if no significant 
dependence on the lattice definition of gauge fields and 
measure. 

Another unambiguous result is the emergence of a 
transverse gluon mass M oc l/a in the strong-coupling 
limit of minimal lattice Landau gauge. Both the gluon 
and ghost propagator show this massive behavior at small 
momenta corresponding to 

Z ~ q 2 /M 2 and G = const. 

for a 2 q 2 <C 1. This massive low-momentum branch of the 
data, however, depends strongly on which lattice defini- 
tion is being used for the gauge fields and their mea- 
sure. This is typical for a mass counter-term on the lat- 
tice and demonstrates the breakdown of lattice Slavnov- 
Taylor identities (STIs) and BRST symmetry in minimal 
lattice Landau gauge beyond perturbation theory. 

It is still possible that this ambiguity disappears in the 
continuum limit, eventually. But because it is a com- 
bination of ultraviolet (mass counter-term) and infrared 
(breakdown of STIs) effects, this might take very fine lat- 
tice spacings in combination with very large volumes and 
therefore who-knows-how big lattices to verify explicitly. 
As we have shown, the ambiguity is definitely present at 
commonly used values of the lattice couplings in SU(2). 

It would obviously be desirable to have a BRST sym- 
metry on the lattice which could then provide lattice 
Slavnov- Taylor identities beyond perturbation theory. 
Non-perturbative lattice BRST has been plagued by the 
Neuberger 0/0 problem, but its improved topological 
understanding provides ways to overcome this problem 
[52l | . It will be particularly interesting to see whether the 
strong-coupling behavior of the propagators will change 
in such approaches and whether this can lead to an unam- 
biguous definition of lattice Landau gauge beyond pertur- 
bation theory and in the strong-coupling limit. Further 
studies and comparisons of different approaches to non- 
perturbative gauge fixing on the lattice in two and three 



dimensions, in particular in the strong-coupling limit will 
thereby be valuable next steps. 



Note added in revised version 

After completion of the original version of this pa- 
per, a follow-up study of the strong-coupling limit ap- 
peared [53[ . There, independent data is presented for 
the SLG gluon and ghost propagators at f3 = in three 
and four dimensions. While the four dimensional data is 
in complete agreement with the corresponding data pre- 
sented here, the discussion of our results in [53| needs 
some clarification. 

First, we have no intention to obscure the fact that 
all simulations of gluon and ghost propagators using cur- 
rent lattice implementations of the Landau gauge provide 
evidence for the qualitative low momentum behavior of 
the decoupling solution (0) in three and four dimensions. 
This includes all currently explored ranges for the lat- 
tice size L and coupling j3, and our (3 = results are no 
exception. 

In fact, in the strong-coupling limit the scaling solution 
would correspond to straight lines in double-logarithmic 
plots as indicated by the dashed/dotted lines in Figs. [TJ 
[5] and 1 101 possibly up to finite- volume corrections which 
would gradually disappear with increasing L/a. This is 
clearly not what we observe. Rather, the plateaus in the 
data at small lattice momenta a 2 q 2 for both, the gluon 
propagator in Fig. [T] and the ghost dressing function in 
Fig. [5l are a clear indication of a massive (decoupling) 
behavior with almost negligible finite-volume effects as 
we have demonstrated. 

It is nevertheless interesting, however, that we do ob- 
serve for the first time on the lattice a scaling behavior Q 
with a scaling exponent k right inside the expected range, 
in the large momentum regime of the strong-coupling 
gluon propagator. This has been confirmed by the data 
of Ref . [53J and more recently also in 42j . The results for 
the ghost propagator are less conclusive but its behavior 
at large lattice momenta tends to approach a form con- 
sistent with the scaling relation (Q}, even though a clear 
signature for scaling over a significant range of momenta 
is not observed. The ghost propagator is most sensitive 
to the treatment of Gribov copies and thus to the nonper- 
turbative completion of the Landau gauge. In fact, there 
is some recent evidence that the Gribov ambiguity can be 
used to tune the ghost propagator to a scaling form in the 
infrared [42|, [H[ . This would identify the one-parameter 
freedom in the continuum solutions to functional equa- 
tions as a second gauge parameter to complete the 
Lorenz condition nonperturbatively. 

It is an unambiguous result that without this tuning 
the strong-coupling data for both propagators is well de- 
scribed by the decoupling form at small momenta in lat- 
tice units, and by scaling at large. To investigate the 
properties of the scaling branch at large lattice momenta, 
we provide three different models to fit this branch with 
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results that are in good agreement of one another, and 
that give a feeling for the systematic fit-model uncer- 
tainties at the same time. In particular, the pure scal- 
ing model must fail when extending the range of the fits 
too far beyond the scaling branch into the small momen- 
tum region, and we see that explicitly happening in the 
corresponding x 2 /ndf in Fig. [3j An alternative is to de- 
fine effective exponents by logarithmic derivatives and to 
monitor their momentum dependence [53] • This requires 
increased statistics but does not change the general the 
properties of the scaling branch and our conclusions [42j ■ 
These does not imply, however, that we suggest to dis- 
miss the small momentum data, as we believe we were 
being misinterpreted in Ref. (53|. 

This is in contrast to the authors of Ref. [53| who write 
in their revised version that the data in the scaling region 
should be discarded. They argue that the scale aM ~ 1 
for the transition between decoupling and scaling should 
be essentially given by gAqcd and that momenta p in the 
scaling region therefore fail to satisfy p <C Aqcd- Even if 
the argument were correct, which it is not as we will ex- 
plain, we would not understand why any data should be 
discarded. It would seem even more surprising and im- 
portant to understand, if the strong-coupling limit were 
to show infrared scaling for momenta on the order of 
Aqcd- Luckily that is not the case, however. Here is the 
flaw in the argument: The authors of [53[ use a recent 
lattice study of the phase diagram of QCD with one fla- 
vor of staggered quarks in the strong-copuling limit 55] 



to assign 'physical' lengths to the lattice spacing which 
correspond to values of oAqcd between 0.62 and 1.92 de- 
pending on which hadron mass or decay constant of the 
one-flavor strong-coupling model is being used to set the 
absolute scale. There is no absolute scale in the strong- 
coupling limit, however, not in the chiral limit of the 
one-flavor model, let alone in the pure gauge theory. For 
the latter, the relevant masses are those of the glueballs 
which all behave as am cx — ln/3 — > oo for j3 —> 0. It does 
make sense to compare mass ratios, but the relation be- 
tween the absolute scale and Aqcd of the scaling region is 
lost at P = 0. In units of the latter, the only correct inter- 
pretation of the strong coupling limit is a — > oo. The fact 
that the strong-coupling mesons and baryons are point- 
like (in lattice units) actually demonstrates that also for 
the quark model. 

We moreover emphasize once more that the massive 
decoupling branch at small lattice momenta is discretiza- 
tion dependent. We believe that this is important and 
needs to be understood. It can hardly be due to vio- 
lations of rotational invariance as investigated, and not 
surprisingly found to be very small, in 53]. Note that we 
did apply the usual cylinder cuts to our data in order to 
reduce those effects as far as possible. Moreover, we only 
considered lattice momenta with a 2 q 2 < 14 in all our fits 
to the scaling branch as a further precaution. To us the 
observed discretization dependences signal a breakdown 
of Slavnov- Taylor identities in minimal lattice Landau 
gauge beyond perturbation theory as we explain. Maybe 
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FIG. 12: The gluon propagator in lattice units versus a 2 q 2 
comparing data for /3 = and j3 — 0.1. The /3 = data are 
the same as in Fig. [T](32 4 lattice). 



we are over-interpreting the strong-coupling limit which 
unveiled this problem, but the hope that it simply disap- 
pears in the scaling region seems overly optimistic. 

The scaling properties such as exponent and coupling, 
on the other hand, appear to be robust under variations 
of the discretization of the gauge fields as far as we can 
judge from the available data. It is also clear, however, 
that the scaling branch gets suppressed as soon as the 
gauge action is turned back on again, as demonstrated in 
Fig. [T2J and that one is left with the decoupling branch 
for low momenta at finite /3. This emphasizes the im- 
portance of understanding any discretization ambiguity 
of the associated gluon mass, before concluding that this 
mass is now firmly established. 
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Appendix A: Standard lattice Landau gauge (SLG) 

In this Appendix we summarize the standard approach 
that is nowadays widely adopted to study the gluon 
and ghost propagators of an SU(2) [or SU(3)] Landau 
gauge theory using lattice simulations. Typically, such 
simulations first generate a well-thermalized ensemble of 
gauge-invariant link configurations which then, in a sec- 
ond step, is fixed to what above we have called the stan- 
dard lattice Landau gauge (SLG). This is most conve- 
niently done by employing either a steepest-descent al- 
gorithm, like, e.g., the Fourier-accelerated gauge-fixing 
of Refs. [5f| |57|, or an over-relaxation algorithm [58| . 
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Those iterative algorithms are designed to relatively fast 
find (for a fixed U) a set of local gauge transformations 
U XIX -> C/| M = g x U xll gl +ji that minimizes the gauge func- 
tional of SLG which for the SU(2) gauge group is of the 
form 



V U [g}=4'£( 



I -itrl/^ 



(Al) 



Any minimum of V;y[<?] automatically ensures that the 
(lattice) Landau gauge condition, i.e., the lattice (back- 
ward) derivative 



F X {A 9 ) 







(A2) 



is satisfied if the gauge-fixed (lattice) gluon fields are 
given through the standard definition 



1 

2ia 



(A3) 



in terms of the gauge-transformed link U!* . Practically, 
in order to fulfill condition (| A2|) with sufficient accuracy, 
the gauge- fixing algorithm is iterated until, e.g., the stop- 
ping criterion 



maxtr[(V^)(V^ + )] < 10" 
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(A4) 



is satisfied for all lattice sites. 

For the sake of completeness we mention that gauge 
functionals, like (IA1[) . do have a huge number of differ- 
ent local minima. The corresponding gauge-fixed config- 
urations do all satisfy Eq. (|A2[) and are related through 
local gauge transformations to each other. This ambi- 
guity, known as the Gribov-copy problem on the lattice, 
causes minimization algorithms to find different gauge- 
fixed configurations all with a different value Vjj [g] . This 
is, in particular, inevitable in the strong-coupling limit. 
The ghost propagator is known to be affected by the 
Gribov-copy problem while there seems to be no influ- 
ence on the gluon prop agator, see, e.g., Refs. [4!l[5(| and 
for (3 = Ref. [U, in particular. 

To calculate the gluon propagator in momentum space, 
Fast Fourier transformations (FFTs) are applied trans- 
forming the gauge-fixed gluon fields A xfl = A® fl a a /2 (the 
superscript g is dropped in what follows) into momentum 
space from which the gluon propagator is then obtained 
for any momenta € (N^/2, N^/2) as the Monte-Carlo 
average 



DfM = (A a Jk)A b v (-k)) l 



(A5) 



over gauge-fixed configurations U. Since for the standard 
Wilson gauge action the propagator's tree-level form on 
the lattice is of the form 



q 2 (k) J q 2 (k) 



1 



where 



9m 0) 



■ sin 



(A6) 



(A7) 



it is natural to associate with the physical momentum. 
Note that then not only the continuum tensor structure 
of the gluon propagator is retrieved but also the gauge 
condition is manifest in momentum space, i.e., 



q M {k)A«{k) = 



(A8) 



The gluon dressing function is straightforwardly obtained 
from the product Z = q 2 D{k) where D(k) = D™(k)/9 
(summing over a and /i). 

Compared to the gluon propagator, determinations of 
the ghost propagator 

c a \k) = £y*(*-io ({M- l )f y ) u 

X y 

are significantly more compute intensive as they involve 
inversions of the Faddeev-Popov (F-P) operator M. On 
the lattice, M is the Hessian of the gauge functional Vu[g] 
and for the SU(2) SLG of the form 



Mfy[U Xtl ]=j2 



+ J s ab s xy 



{<J ah + £ abj <j6 x+ ^ y (A9) 



where we have used the same notation as in Eq. (|12|) and 
assumed periodic boundary conditions. M for the SU(3) 
SLG can be found, e.g., in [Hof . 

Due to its zero eigenvalues, with corresponding con- 
stant eigenmodes, it is impossible, however, to construct 
the inverse of M . Nevertheless, all is not lost, as for 
the final analysis one is usually interested only in data of 
G at finite momenta for which M can be inverted on a 
(non-constant) vector of plane waves £° (c, k) — S ac e 27Tlkx 
with k ^ 0. Typically, the momenta k are chosen to 
survive both the cylinder and cone cut 60] and then 
M X yCy(c,k) = Q{c,k) is solved for ((c,k) yielding the 
ghost propagator G ab (k) = S ab G(k) (directly in momen- 
tum space) where G(k) is the MC average 



1 3 

G(k) = -J2m,c)-ak,c)) l 



(A10) 



c=l 



In this way, also translational invariance is used to its full 
capacity reducing the statistical noise to a minimum. 

For the inversion it is highly recommended to use a 
pre-conjugate gradient algorithm, e.g., that of Ref. [50l |. 
This precondition has been proven to drastically accel- 
erate computations, in particular, when lattices sizes are 
as big as we have used for this study and can also be 
straightforwardl y appli ed to the case of Coulomb gauge 
as done, e.g., in [6lLl62|. For further details on this tech- 
nique refer to [631 ] . 



15 



At tree-level the F-P operator is simply — 5 ab A xy where 
A denotes the lattice Laplacian in four dimensions. Con- 
sequently, the tree-level form of the ghost propagator on 
the lattice is 



G a b (k) 



where q M is that of Eq. (|A7|> . 



:a b 



1 



q 2 (k) 



(All) 



Appendix B: Modified lattice Landau gauge (MLG) 

In this Appendix we summarize the Modified lattice 
Landau gauge (MLG) of Ref. [43[ and provide the reader 
with the necessary information for the gauge group 
SU(2) as used in Sec. El 

The MLG is a novel implementation of the Landau 
gauge for lattice gauge theories which was introduced in 
Ref. [43[ for the gauge groups U(l) and SU(2). It is based 
on stereographic projection to define lattice gauge fields 
and extensions to the gauge group SU(3) or to Coulomb 
gauge are also possible. 

When comparing MLG to the ever popular SLG, there 
is no advantage that the SLG has over the MLG. A 
promising particular feature of the MLG on the other 
hand is that it provides a way to perform gauge- fixed MC 
simulations sampling all Gribov copies of either sign (of 
the Faddeev-Popov determinant) in the spirit of BRST. 
This is because in MLG there is no perfect cancellation of 
Gribov copies of opposite sign, known as the Neuberger 
0/0 problem, which prevented us from performing such 
simulations for 20 odd years. Note that in the standard 
lattice approach to Landau gauge, i.e., via a minimiza- 
tion of the SLG functional [Eq. (|A1[) ] . only Gribov copies 
within the first Gribov region, i.e., with positive sign, are 
sampled. Gauge-fixed Monte Carlo simulations via MLG, 
however, would enable us to sample beyond that region. 
Indeed such simulations are not meant to compete in any 
sense with standard MC simulation of lattice gauge the- 
ory, as gauge-invariant observables would be completely 
unaffected by that, but rather to provide a theoretical 
sound framework for studying non-perturbative proper- 
ties of gauge- variant correlation functions. 

This particular feature of MLG will be explored in a 
forthcoming study. Here we simply have used the MLG 
for comparison in the standard way, i.e., we gauge- fix 
configurations via minimization of the MLG functional, 



Vu\g] 



-8 



5>G 



(Bl) 



Consequently, only the first Gribov region of MLG is 
sampled, however when doing so we are in the fortu- 
nate position of having two alternative lattice implemen- 
tations of SU(2) Landau gauge theory which differ at 
finite lattice spacing, a, but meet in the continuum limit. 
When comparing SLG to MLG data the impact of dis- 
cretization errors can then be seen at any finite a. 

Leaving aside the issue of Gribov copies, we adapt 
the Fourier- acceleration of Ref. (5(| to gauge-transform 
configurations such that they satisfy the lattice Landau 
gauge condition, here that of MLG, 







(B2) 



This automatically maximizes the MLG functional 
[Eq. (1BT]) ] when A xfl , the lattice gluon field of MLG, is 
given through 



where 



E4 M 



1 + 5 tr U XI1 



(B3) 



(B4) 



In our study the Fourier-accelerated gauge-fixing is iter- 
ated until the stopping criterion 



max tr 



V h „A x ^ b ,Ai 



Xjl 



< 10 



-13 



(B5) 



is met at all lattice sites. 

Note that gauge-fixed configurations satisfying 
Eq. flB5j| do not satisfy that of SLG [Eq. (lA4|) ] and vice 
versa. Nevertheless, transversality of the corresponding 
lattice gluon field is ensured in both cases, if momenta 
are associated with Eq. (|A7|) and the standard midpoint 
definition is assumed. That is, 



£>(A)A»(*) = 



(B6) 



The gluon propagator of MLG is straightforwardly 
constructed as in SLG with A substituted through A. 
Similar holds for the ghost propagator, though in MLG 
the F-P matrix is of the form 



X.fl 



16 



Mab \ * ) ( ~0 rafc ■ abc~c , ^ -a -b \ c 

$xy 



1 2 x—ji^^x—jXtfJ. 



X — fl,fl 



ab _ ^abcjjC 



+ k 



&X — LL. 



(B7) 



where tt° = 2u°/(l + u°) and u = 2u/(l + u°). Due to 
the logarithm in Eq. ([Bill there are M a M f> -terms quadratic 
in the projected variables u. Apart from those, M is of 
the same form as the standard F-P operator in SU(2) 



[Eq. (|A9[) ] with the u's replaced by u. Therefore, the tree- 
level forms of the gluon and ghost propagators, Eq. (|A6[) 
and (|A11|) . are also valid in MLG. 
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